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Abstract. Using the methods of quantum trajectories we study effects of dissipative decoherence on the 
accuracy of the Grover quantum search algorithm. The dependence on the number of qubits and dissipation 
rate are determined and tested numerically with up to 16 qubits. As a result, our numerical and analytical 
studies give the universal law for decay of fidelity and probability of searched state which are induced by 
dissipative decoherence effects. This law is in agreement with the results obtained previously for quantum 
chaos algorithms. 
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Nowadays the quantum computing attracts a great in- 
terest of the scientific community pQ. The main reason of 
that is due to the fact that certain quantum algorithms 
allow to perform computations much faster than the usual 
classical algorithms. The famous example is the Shor algo- 
rithm which performs factorization of integers on a quan- 
tum computer exponentially faster than any known classi- 
cal factorization algorithm 2 . However, at present there 
is no mathematical prove about efficiency of potentially 
possible classical algorithms that in this case gives cer- 
tain restrictions on the comparative efficiency of classical 
and quantum algorithms. The situation is different in the 
case of the Grover quantum search algorithm |3] . Indeed, 
it has been proved that it is quadratically faster than any 
classical algorithm (see e.g. £Q and Refs. therein). 

In addition to the question of quantum algorithm effi- 
ciency it is also important to know what is the accuracy 
of a quantum algorithm in presence of realistic errors and 
imperfections. The accuracy can be characterized by the 
fidelity / |2] defined as a scalar product of the wave func- 
tion of an ideal quantum algorithm and the wave func- 
tion given by a realistic algorithm (see e.g. (T|). In gen- 
eral it is possible to distinguish three types (classes) of 
errors. The first class can be viewed as random unitary 
errors in rotational angles of quantum gates. This is the 
mostly studied case which has been also analyzed for var- 
ious quantum algorithms with the help of numerical sim- 
ulations with up to 28 qubits (see e.g. 5 6,7,8,9,10 ). It 
has been shown that in such a case the fidelity decays 
exponentially with the number of quantum gates n g and 
with a rate 7 which is proportional to a mean square of 
fluctuations in gate rotations. The second class of errors is 
related to static imperfections (static in time). They are 
produced by static residual couplings between qubits and 



static energy shifts of individual qubits which may gen- 
erate many-body quantum chaos in a quantum computer 
hardware jllj . This second type of errors, e.g. static im- 
perfections, gives a more rapid decay of fidelity as it has 
been shown in |12Ul()j . In the case of quantum algorithms 
for a complex dynamic these imperfections lead to the fi- 
delity decrease described by a universal decay law given 
by the random matrix theory jl(J|. The two former classes 
are related to unitary errors. However, there is also the 
third class which corresponds to the case of nonunitary 
errors typical to the case of dissipative decoherence. This 
type of errors has been studied recently for the quantum 
baker map jlHj and the quantum sawtooth map |14| . It 
has been shown that the exponential decay rate of fidelity 
is proportional to the number of qubits. The dissipative 
decoherence is treated in the frame of Lindblad equation 
for the density matrix ^5]. A relatively large number of 
qubits can be reached by using the methods of quantum 

The quantum algorithms studied in 7 8 Till 1011 12lll.il 

describe quantum and classical evolution of dynami- 
cal systems. However, it is also important to analyze the 
accuracy of realistic quantum computations for more stan- 
dard algorithms, e.g. for the Grover algorithm. The effects 
of unitary errors of the first and second classes have been 
studied in 123] and |2U respectively. It has been shown 
that the accuracy of computation is qualitatively different 
in case of random errors in gates rotations |23j and in the 
case of static imperfections [21]. Thus it is important to 
analyze the effects of dissipative decoherence in the Grover 
algorithm to have a complete picture for this well known 
algorithm. For that aim we use the approach developed in 
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To study the effects of dissipative decoherence in the 
Grover algorithm we use the same notations as in |24| . 
Thus, the computational basis of a quantum register with 
N = 2 n " states ({\x)}, x = 0, . . . , TV - 1) is used for the 
algorithm itself. According to PE], the initial state |f/ , o) = 
Y^Zq \ X )/VN, is transfered to the state 

fo&(t)) = G%>> m 

= sin ((t + 1/2)w g )|t) + cos ((t + l/2)w G )|r?) , 1 J 

after t applications of the Grover operator G. Here, loq = 
2 arcsin( N) is the Grover frequency, \t) is the search 
state and \rf) = J^i^r^ \ x )/°JN~— 1- Hence, the ideal 
algorithm gives a rotation in the 2D plane (|t), |t^)). One 
iteration of the algorithm is given by the Grover operator 
G and can be implemented in n g = Ylntot ~ 42 elemen- 
tary gates including one-qubit rotations, control-NOT and 
Toffolli gates as described in [21]. The implementation of 
all these gates requires an ancilla qubit so that the total 
number of qubits is n to t — n q + 1. 

To study the effects of dissipative decoherence on the 
accuracy of the Grover algorithm we follow the approach 
used in |13U14| . The evolution of the density operator p(t) 
of open system under weak Markovian noises is given by 
the master equation with Lindblad operators L m (m = 
1, • • • ,n to t) 

P = -\ [HeffP - pH\ f f ]+J2 L mpL\n , (2) 

m 

where the Grover Hamiltonian Hq (G — exp(-iHc)) 
is related to the effective Hamiltonian H e ff = Hq — 
iti/2'Y Jm L\ n L m and m marks the qubit number. In this 
paper we assume that the system is coupled to the en- 
vironment through an amplitude damping channel with 
L m = a m ^/T, where a m is the destruction operator for 
m— th qubit and the dimensionless rate r gives the decay 
rate for each qubit per one quantum gate. The rate r is 
the same for all qubits. 

This evolution of p can be efficiently simulated by av- 
eraging over the M quantum trajectories (see e.g. 
H7llT8 ,19 .20 , 2*ni22| ) which evolve according to the follow- 
ing stochastic differential equation for states \ip a ) (a — 
M): 

-LlL m )\r)dt + £ m lj=Q^ - 1) mdN m , 

(3) 

where (}^ represents an expectation value on \ip a ) and 
dN m are stochastic differential variables defined in the 
same way as in ^5] ( see Eq.(10) there). The above equa- 
tion can be solved numerically by the quantum Monte 
Carlo (MC) methods by letting the state \tp a ) jump to 
one of L m \ijj a ) /\L m \ijj a )\ states with probability 
dp m = \L m \ip a )\ 2 dt [T33 or evolve to 

(1 — iH e f fdt/fr)\il} a ) I ' y/\ — J2 m dp m with probability 1 — 
^ m dp rn . Then, the density matrix can be approximately 
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Fig. 1. (color online) Probability of a searched state wg(£) 
(top), the weight of 4-state subspace un{t) (middle) and the 
fidelity f(t) (bottom) in the Grover algorithm as a function of 
the iteration step t for n tot — 12 qubits, to = 35.5. Dotted 
curve shows data for the ideal Grover algorithm, and solid 
curves from top to bottom correspond to one qubit decoherence 
rate r = (1,2,4,8) x 10~ 5 . Number of quantum trajectories 
used in simulations is M — 400. 

expressed as 

1 M 

P {t) « (\m)W)\) M = jjjEhrww ®! - w 

a=l 

where ()m represents an ensemble average over M quan- 
tum trajectories \ip a (t)}. Hence, an expectation value of 
an operator O is given by (O) — Tr(Op) w (0)m- 

In presence of dissipative decoherence the fidelity / of 
quantum algorithm is defined as 

f(t) = (Mt)\ P (t)\Mt)) « ^El^oWI^W)! 2 . ( 5 ) 

a 

where \^o{t)) is the wave function given by the exact al- 
gorithm and pit) is the density matrix of the quantum 
computer in presence of decoherence, both are taken af- 
ter t map iterations. Here, pit) is expressed approximately 
through the sum over quantum trajectories (see also [Dj]1. 

The effects of dissipative decoherence for the Grover al- 
gorithm with n to t = 12 qubits are presented in Fig.^ The 
data for the probability of the searched state wg it) clearly 
show that the Grover oscillations, which have the period 
2tc ~ tt/^Gj start to decay with the number of iterations 
t due to decoherence. It is not so easy to extract an expo- 
nential decay superimposed with oscillations. Therefore, it 
is convenient to study also the decay of the 4-states prob- 
ability w±it) (| r) and 1 77) for two states of ancilla qubit, 
see |23|). This probability has a pure exponential decay 
with the rate 7: w&it) = exp(— jt). The fidelity decay has 
the same rate 7 (see Fig. 

It is interesting to note that the decay of Grover oscil- 
lations depends on the binary representation of searched 



O.V.Zhirov and D.L.Shepelyansky: Dissipative decoherence in the Grover algorithm 



3 




T 



1 2 

t/t G 

Fig. 2. (color online) Probability of a searched state wa{t) in 
the Grover algorithm as a function of the iteration step t for 
titot = 9 qubits. Solid curves from top to bottom correspond to 
different numbers of units in the binary expansion of the num- 
ber of a searched state, n u = 1,2,4,6,8. Number of quantum 
trajectories used in simulations is M — 2000. 




Fig. 3. Husimi function in the Grover algorithm at different 
number of iterations: t = 1, 17, 35, from left to right. Top panels 
correspond to the ideal Grover algorithm, and bottom panels 
correspond to the single-qubit decay rate F — 2 • 10~ 4 (n to t = 
12). The horizontal axis shows the computational basis x = 
0, 2N — 1, while the vertical axis represents the conjugated 
momentum basis. Density is proportional to grayness changing 
from maximum (white) to zero (black). 



state \t). Indeed, the larger is the number n u of spin-up 
states (numer of units) in the binary representation of |r) 
the faster is the decay of oscillations as it is shown in 
Fig. El The physical origin of this effect is quit clear since 
in the amplitude decoherence channel the decay r takes 
place only for spin-up qubits. However, even if this effect is 
clearly seen in Fig. [2 in average it is not very strong since 
during evolution the qubit rotates between two states that 
leads to averaging of this effect. In the following we will 
neglect the small deviations produced by this effect and 
will analyze the averaged behaviour. 
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Fig. 4. (color online) Effective rate 7 for decay of the prob- 
ability of searched state (top), effective rate 74 for decay of 
the weight of 4-state invariant Grover subspace (middle) and 
effective fidelity decay rate 7/ (bottom) as a function of a sin- 
gle qubit decay rate F; n e ff = n g ntot- Data correspond to 
ntot = 8 (open squares), 9 (open triangles), 10 (open circles) 
and to ntot = 12 (full squares), 15 (full triangles), 16 (full 
circles). Full lines show the average dependence @. 

A pictorial image of the algorithm accuracy can be ob- 
tained with a help of the Husimi function 25,24 which 
is shown in Fig. EI ln the ideal algorithm the total proba- 
bility is transfered from the initial state (horizontal white 
half line corresponding to one state of ancilla qubit) to 
the searched state (vertical white line) after 35 iterations. 
In presence of dissipative decoherence the probability of 
searched state is significantly reduced and probability is 
transfered to the state \rj) and many other states of the 
computational basis (Fig. EJ) . It is interesting to note that 
these transitions give homogeneous distribution in the com- 
putational basis (at one state of ancilla qubit) and select 
specific momentum states in the momentum representa- 
tion that mainly corresponds to dissipative flips of qubits 
to zero state. 

From the data similar to those of Fig. ^ it is pos- 
sible to obtain the exponential decay law for WG,iVi,f 
(oc exp(— jt)) and find from it the corresponding decay 
rates 7,74,7/. Their dependence on the one qubit decay 
rate r is shown in Fig. El for n to t — 8, ....16. In all three 
cases the dependence is well described by the relation 

7 = Crn g n t ot , (6) 

where C ~ 0.4 is a numerical constant. This result gives 
the decay rate 7 which is close to the maximal decay rate 
n to t which corresponds to the state with all 
qubits in the up state. Quite naturally 7 is proportional 
to the total number of qubits ntot since a flip of each qubit 
from up state gives the decay of fidelity for the whole wave 
function. The constant C is not sensitive to a variation of 
the number of qubits that is illustrated in Fig. El Indeed, 
a fit for two groups of qubits [n to t — 8, 9, 10 and n to t — 
12, 15, 16) gives close values of C (see Fig. El- We note that 
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Fig. 5. (color online) Ratios of effective decay rates 7, 74 and 
7/ (as in Fig|lJ to the maximum decay rate r m ax = rn g ntot 
as a function of r. The dashed and full lines give the fit values 
of constant C obtained for two groups of qubits (dashed for 
ntot = 8, 9, 10 and full for ntot = 12, 15, 16). The symbols are 
the same as in Fig0] 



C is close to the value 1/2 since in average only a half of 
qubits is in the up state that reduces the value of 7 by a 
factor 2 comparing to the maximum decay rate r max . 

The result JTjJ) is in agreement with the dependence 
found in |13U14| for other quantum algorithms with dis- 
sipative decoherence. This means that the decay rate re- 
lation (JBJ gives a universal description of dissipative de- 
coherence in various quantum algorithms. Therefore it is 
possible to compare the three classes of quantum errors 
described at the beginning. The comparison shows that 
the most rapid decrease of fidelity, and thus the accuracy 
of quantum computation, is produced by static imperfec- 
tions. Thus it is necessary to develop specific error correc- 
tion methods which will be able to handle effects of static 
imperfections in quantum algorithms. First steps in this 
direction are done recently in |26U27U28| . 

This work was supported in part by the EC IST-FET 
projects EDIQIP and EuroSQIP and (for OVZ) by RAS 
Joint scientific program "Nonlinear dynamics and soli- 
tons" . 
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